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We present a convergent variational basis-set calculational scheme for elas- 

. tic scattering of positronium atom by hydrogen atom in S wave. Highly cor- 
\D _ 

CO ' related trial functions with appropriate symmetry are needed for achieving 
O ' 

cn ■ convergence. We report convergent results for scattering lengths in atomic 

units for both singlet (= 3.49 ± 0.20) and triplet (= 2.46 ± 0.10) states. 
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Lately, there has been interest in the experimental |jj and theoretical |2|-§| studies of 
ortho positronium (Ps) atom scattering by different neutral atomic and molecular targets. 
The Ps-H system is theoretically the most simple and fundamental and a complete under- 
standing of this system is necessary before a venture to more complex targets [|7|-|9|]. There 
have been R-matrix |J, close-coupling (CC) HJlO| and model-potential calculations 
for Ps-H scattering. Here we present a convergent variational basis-set calculational scheme 
for low-energy Ps-H scattering in S wave below the lowest Ps-excitation threshold at 5.1 eV. 
Using this method we report numerical results for scattering length of electronic singlet and 
triplet states. 

A recent study based on a regularized nonlocal electron-exchange model potential 
yielded low-energy (total) cross sections in agreement with experiment for Ps scattering 
by He ||, Ne |7|], Ar [|7| and H 2 ||. For the Ps-H system the model-potential results for 
S-wave singlet binding and resonance energies are in agreement with accurate variational 



estimates |L2] . It would be interesting to see if the model-potential result for the Ps-H singlet 
scattering length agrees with the present convergent calculation. 

Because of the existence of two identical fermions (electrons) in the Ps-H system, one 
needs to antisymmetrize the full wave function before attempting a solution of the scattering 
problem. The position vectors of the electrons — ri (Ps) and r2 (H) — and positron (x) with 
respect to (w.r.t.) the massive proton at the origin are as shown in figure 1. We also use the 
position vectors Sj = (x + r,)/2, pj = x — r,-, j = 1, 2, r 12 = — r 2 . The fully antisymmetric 
state ip£ of Ps-H scattering is given by = -Ail^k) = (1 i P^l^k) = l^k) =t l^k) where 
k is the incident momentum, the antisymmetrizer A\ is (1 + Pi 2 ) for the singlet state and 
(1 — P12) for the triplet state with P 12 the permutation operator of electrons 1 and 2. The 
function i/j^ refers to the Ps-H wave function with electron 1 forming the Ps as in figure 1 
and ipl. refers to the same with the two electrons interchanged. 

The full Ps-H Hamiltonian H can be broken in two convenient forms as follows H = 
Hi + V\ = H 2 + V2 where Hi includes the full kinetic energy and intracluster interaction 
of H and Ps for the arrangement shown in figure 1 and Vy is the sum of the intercluster 
interaction between H and Ps in the same configuration, if 2 and V 2 refer to the same 
quantities with the two electrons interchanged: 



Vi 



1 1 1 

x ri r l2 



1 

Pi 



V 2 = P 12 V, 



1 1 1 

x r 2 r l2 



1 

Pi 



The fully antisymmetric state satisfies the Lippmann-Schwinger equation [13 



(2) 



where the channel Green's function Gi = (E + iO — Hi)" 1 and the incident wave [0^) satisfies 
(E — Hi)\(f)\) = 0. The incident Ps energy = 6.8k 2 eV. We are using atomic units (au) in 
which ao = e = m = h = l, where e (m) is the electronic charge (mass) and ao the Bohr 
radius. Using the definition of the antisymmetrized state we rewrite (Q) as [ff3| 



(3) 



M a = Ki^i + (£ - #i)(l - A) = A^i + (1 - - Hi). 
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The properly symmetrized transition matrix for elastic scattering is defined by \ T A | = 

(^kl^ilV'k') = (^kl^i-^-ilV'k) = (V'kl^i^il^k) EH- ^ basis-set calculational scheme for the 
transition matrix can be obtained from the following expression [IT? 



M T A \<fc) = (i&lAiVM) + (KlAiVM) ~ MAM - MiGiAiVi\ipl). 
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Using (|3|), it can be verified that @ is an identity if exact scattering wave fumctions 
are used. If approximate wave functions are used, (|5]) is stationary w.r.t. small variations 
of \ifj^) but not with (V'kl- This one-sided variational property emerges because of the lack 
of symmetry of the formulation in the presence of explicit antisymmetrization operator A\. 
However, this variational property can be used to formulate a basis-set calculational scheme 
with the following trial functions |T3 
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Substituting (§) into @ and using this variational property w.r.t. [■?/>£) we obtain Jl4 

N 
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(fmlAxVt - [AiVx + (1 - Ai)(E - HJlGtAiVxlfn). 



(8) 



Using the variational form (|7|) and definition ((f^T \<p^) = (V'kl-^-i^-l^k) we obtain the 
following basis-set calculational scheme for the transition matrix 

N 

«\T A \<Pl)= E (tilAiVilMD^iUlAMlti). (9) 

m,n=l 

(§) and (|9]) are also valid for the K matrix and in partial-wave form where the momentum- 
space integration over the Green's function Gi should be performed with the principal- value 
prescription. 

In the calculation, the basis functions are taken in the following form 

/m(r 2 ,Pi,Si) = ^(r 2 )r/(p 1 )e^ T2 - a - pl -^ Sl - 7 - (p2+ri2) -^ (s+ri) Sm(A:Sl) (10) 

ksi 

where (p(r) = exp(— r)/y/n and r](p) = exp(— 0.5p)/v87r represent the H(ls) and Ps(ls) 
wave functions, respectively. For elastic scattering the direct Born amplitude is zero and the 
exchange correlation dominates scattering. To be consistent with this, the direct terms in the 
form factors (/ m |^i^i|0k) an d (0kl^i^i|/n) are zero with the above choice of correlations in 
the basis functions via 7 m and p m . This property follows as the above function is invariant 
w.r.t. the interchange of x and ri whereas the remaining part of the integrand in the direct 
terms changes sign under this transformation. A proper choice of the correlation parameters 
7 m and p m is crucial for obtaining good convergence. 

In the following we specialize to the K-matrix formulation in S wave at zero energy, when 
sin(ksi)/(ksi) = 1 in fllPl) . The useful matrix elements of the present approach are explicitly 
written as [|14|] 

| p Sill i pS'2 l 

{4> x p \AxVx\f n ) = ±— / (p(ri)r)(p 2 ) [K]/n(r 2 ,Pi,Si)dr 2 dpidsi (11) 

Z7T J pS2 

1 f sin j)S\ 

(fmlAxVM) = ±— / / m (ri,p 2 ,s 2 )[Vi]v?(r 2 )?7(pi) dr 2 dp 1 ds 1 (12) 

y lix J psi 

(fm\AiVi\f n ) = ±— J f m (ri, p 2 ,s 2 )[Vi)f n (r2, px,si)dr 2 dpidsi (13) 
(fmlMidAiVilfn) w — / dpiUlAiVMlHfilAMlU) (14) 
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where the so called off-shell term (1 — A\){E — Hi) has been neglected for numerical sim- 
plification in this calculation. This term is expected to contribute to refinement over the 
present calculation. In this convention the on-shell K-matrix element at zero energy is the 
scattering length: a = ((pl\K A \(pl) . 



All the matrix elements above can be evaluated by a method presented in [15|. We 
describe it in the following for ((j)p\AiVi\f n ) of ([TT|). By a transformation of variables from 
(r 2 , Pi,&i) to (si, s 2 , x) with Jacobian 2 6 and separating the radial and angular integrations, 
the form factor (]TT|) is given by 
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e - { ar 1+ b Pl /2) e -(cr 2 +d P2 /2) e - lnri2 [V^^df (15) 



where a = 1 + p n , b = 2a n + 1, c = 1 + 6 n and d = 2^ n + 1. Recalling that = 2sj — x, 
ri2 = 2(si — S2), Pj = 2(x — Sj),j = 1,2, we employ the following expansions of the 
exponentials in (O) 
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e -a|2s-x|-6|x-s| = - J2G\ a ' b \s,x)YC m (s)Y lm (x) 
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where the 3^ m 's are the usual spherical harmonics. Using (|T6|) — (p0|) in (|15|) we get 

(0jiAv 1 |/ n > = ±2 8 r e -^d 8l rds 2 s -^^- r dxe-^^i+i) 
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where the Z-sum is truncated to I = L. This evaluation avoids complicated angular integra- 
tions involving Si, S2 and x. These integrals take a simple form requiring straightforward 
numerical computation of certain radial integrals only, which must, however, be carried out 
carefully. The functions Gi, Ji, Ki etc. are easily calculated using fll6|) — (|2( 



G<f> b \s,x) 



(22) 



where Pi(t) is the usual Legendre polynomial and t is the cosine of the angle between s and 
x. The integrals ([12|) and (|13|) can be evaluated similarly. For example 
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where e = 1 + 5 m + fi n , f = 2a n + 2^ m + 1, g = 1 + 6 n + fi m , h = 2a m + 2^ n + 1 and 

We tested the convergence of the integrals by varying the number of integration points in 
the x, Si and s 2 integrals in fl2~T|), (p3|) and ( p4j) and the t integral in (0). The evaluation of 
(pl|), ( p3| ) and (23) essentially involves four-dimensional integration, which is performed with 
caution. The x integration was relatively easy and 20 Gauss- Legendre quadrature points 
appropriately distributed between and 16 were enough for convergence. In the evaluation 
of integrals of type (^) 40 Gauss- Legendre quadrature points were sufficient for adequate 
convergence. The convergence in the numerical integration over S\ and s 2 was achieved 
with 300 Gauss- Legendre quadrature points between and 12. The maximum value of I 



in the sum in fl2~l"D , (^) and (|24]) , L, is taken to be 6 which is sufficient for obtaining the 
convergence with the partial- wave expansions ([16|) — (|20|). 



Table 1: Singlet (a s ) and triplet (at) Ps-H scattering lengths for different L and N. 
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a t a s 
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4.00 3.61 


3.66 4.00 


3.07 3.75 


2.92 3.74 


7 


3.98 4.11 


3.25 3.88 


2.99 4.06 


2.85 4.00 


8 


3.93 4.12 


3.35 3.96 


2.66 3.72 


2.54 3.73 
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3.97 4.15 


3.42 3.83 


2.65 3.92 


2.55 4.06 


10 


3.97 4.22 


3.42 3.92 


2.57 4.42 


2.50 3.45 


11 


4.01 3.82 


3.44 3.91 


3.55 3.75 


2.67 3.80 


12 






2.54 3.72 


2.46 3.73 


13 






2.48 3.47 


2.46 3.49 



In a numerical calculation a judicial choice of the parameters in (|10"D is needed 
for rapid convergence. As the variational method does not provide a bound on the 
result, the method could converge to a wrong scattering length if an inappropriate 
(incomplete) basis set is choosen. After some experimentation we find that to ob- 
tain proper convergence the parameters 5 n and a n should be taken to have both 
positive and negative values, 7„ and fi n should include values close to unity and 
P n should have progressively increasing values till about 1.5. The results reported 
in this work are obtained with the following parameters for the functions f n ,n = 
1, 13: {S n ,a n ,p n ,j n ,fx n } = {-0.5,-0.25,0.3,0.01,0.02}, {-0.5,-0.25,0.5,0.04,0.02}, 
{-0.5, -0.25, 0.7, 0.03, 0.06}, {-0.2, -0.1, 0.6, 0.2, 0.2}, {-0.1, 0.1, 0.8, 0.25, 0.25}, 
{0.2, -0.2, 0.6, 0.35, 0.35}, {-0.1, -0.1, 0.7, 0.4, 0.4}, {0.15, 0.2, 0.8, 0.5, 0.5}, 
{0.12,-0.12,1,0.7,0.7}, {0.2,0.2,1.2,0.9,0.9}, {0.1,0.2,1.3,1,0.7}, {0.2,0.1,1.4,0.7,1}, 
{0.3,0.15,1.5,1,1} 

In table 1 we show the convergence pattern of the present calculation w.r.t. the number 
of partial waves L and basis functions iV used in the calculation. The convergence is smooth 
with increasing L. However, as the present calculation does not produce a bound on the 



result, the convergence is not monotonic with increasing N. The lack of an upper bound in 
this calculation is clearly revealed in table 1 where the results do not decrease monotonically 
as the number of terms in the trial wave function is increased. The unbounded nature of 
the results is consistent with the noted oscillation of the scattering lengths as N increases. 
The oscillation is larger in the singlet state where it is more difficult to obtain convergence. 
Although, the final result for the largest N and L is supposed to be the most accurate, it 
is not known whether this result is larger or smaller than the exact one. Also, an estimate 
of error of this result is not known. It is difficult to provide a quantitative measure of 
convergence. However, from the noted fluctuation of the results for large N for L = 6 
we believe the error in the triplet scattering length to be less than 0.10 and in the singlet 
scattering length to be less than 0.20. The final results of the present calculation are those for 
N = 13 and L = 6 with the above estimate of error: a s = 3.49 ±0.20 au and 04 = 2.46 ±0.10 
au. 

The maximum number of functions (N = 13) used in this calculation is also pretty 
small, compared to those used in different Kohn-type variational calculations for electron- 
hydrogen (N = 56) |TE| and positron-hydrogen (N < 286) |T7| scattering. Because of the 
explicit appearance of the Green's function, the present basis-set approach is similar to the 
Schwinger variational method. Using the Schwinger method, convergent results for electron- 
hydrogen |TB[ and positron-hydrogen [HJ scattering have been obtained with a relatively 
small basis set (N ~ 10). These suggest a more rapid convergence in these problems with a 
Schwinger-type method. 

Now we compare the present results with those of other calculations. While a static- 
exchange calculation by Hara and Fraser [|K| yielded a s = 7.28 au and a t = 2.48 au a 
model potential calculation by Drachman and Houston |T]J produced a s = 5.33 au and 
a t = 2.36 au. A 22-state R-matrix calculation by Campbell et al. [Q] and a six-state CC 
calculation by Sinha et al. j|] yielded a s = 5.20 au, a t = 2.45 au and a s = 5.90 au, a t = 2.32 
au, respectively. These triplet scattering lengths are in qualitative agreement with the 
present result: at = 2.46 ±0.10. However, the singlet scattering lengths of these calculations 
have not yet converged. The disagreement in the singlet channel shows that it is more 
difficult to get the converged result in the attractive singlet channel than in the repulsive 



triplet channel. This is consistent with the common wisdom that a scattering model at low 
energies is more sensitive to the detail of an effective attractive interaction than to that of a 
repulsive interaction. The nonconvergence of these results for the singlet scattering length 
was conjectured before ||. Using a correlation between the S-wave scattering length and 
binding energy for the Ps-H system, the value a s = 3.5 au was predicted in that study || 
in excellent agreement with the present result: a s = 3.49 ± 0.20. 

To summarize, we have formulated a convergent basis-set calculational scheme for S-wave 
Ps-H elastic scattering below the lowest inelastic threshold using a variational expression for 
the transition matrix. We illustrate the method numerically by calculating the singlet and 
triplet scattering lengths: a s = 3.49 ± 0.20 au and a t = 2.46 ± 0.10 au. 
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Figure Caption: 

1. Different position vectors for the Ps-H system w.r.t the proton (p) at the origin in 
arrangement 1 with electron 1 forming the Ps. 




Figure 1 



